Biophysical Journal

Dynamic conformational changes in intracellular domains fundamentally affect the function and pharmacology of many membrane proteins. For example, sodium influx through the cardiac voltage-gated sodium channel (NaV1.5) is rapidly terminated through conformational changes that result in pore closure, a transition known as inactivation. Inactivation involves NaV1.5 intracellular regions, particularly the DIII-DIV linker containing the IFM particle (Isoleucine-Phenylalanine-Methionine) and its dysfunction is a major cause of cardiac arrythmias. However, the conformational changes involved in inactivation and their modulation by auxiliary proteins and clinically used drugs remain incompletely characterized, partly because live-cell, site-specific labeling of intracellular regions with small fluorescent dyes remains challenging. Here, we combine live-cell protein semi-synthesis with voltage-clamp fluorometry (VCF), to track intracellular conformational dynamics of the cardiac sodium channel NaV1.5 and monitor their voltage dependence and kinetics in real time. We identify intracellular conformational changes involved in both fast and steady-state inactivation of NaV1.5 and show that both lidocaine and auxiliary proteins affect the kinetics of conformational changes of the DIII-DIV linker. Our work establishes the combination of protein semi-synthesis and voltage-clamp fluorometry as a powerful approach to dissect intracellular conformational changes in membrane proteins.

Diffusion within the cytoplasm is fundamental to numerous biological processes. Fluorescence recovery after photobleaching (FRAP) is one of the most common method for quantifying molecular diffusivity in living cells using standard laser scanning confocal microscopy (LSCM). However, accurately measuring fast cytoplasmic diffusion (typically >10 m^2/s) is challenging due to rapid recovery kinetics, weak signal-to-noise ratios, post-bleach signal artifacts, and spatial restrictions affecting normalization. While individual challenges have been addressed in specific contexts, a simple and robust framework to quantify cytoplasmic diffusivity remains elusive. Here, we present a FRAP methodology specifically designed to overcome these obstacles. By utilizing the Gaussian function -- the impulse response (ImpRes) of the diffusion equation in an infinite medium -- our approach leverages the full spatiotemporal dataset through a single-equation three-parameter fitting procedure, thus releasing restrictions to small regions of interest and arbitrary initial time-points. The methodology was validated on three datasets of increasing complexity: in silico simulated recovery profiles, in vitro data from FITC-dextran in glycerol solution, and live-cell imaging of free cytoplasmic GFP. Systematic comparison with existing models demonstrates that the ImpRes approach significantly reduces sensitivity to noise and imperfect fluorescence normalization, while remaining robust against short-term biases, such as transient probe photo-activation. Given its robustness under realistic experimental conditions and its ease of implementation, the proposed FRAP methodology provides a reliable tool for quantitative cytoplasmic analysis.

During the influenza viral life cycle, the viral glycoprotein hemagglutinin (HA) mediates binding, entry, and fusion. Densely packed clusters of HA trimers at the plasma membrane are required to produce infectious virions; however, the mechanism of HA clustering is still unknown. We have shown previously that HA co-clusters with and modulates phosphatidylinositol 4,5-bisphosphate (PIP2) in host cell plasma membranes (PM). Here, we further characterize the relationship between HA and PIP2 using molecular dynamics simulations (MD) and fluorescence photoactivation localization microscopy (FPALM) to elucidate a mechanism of HA-PIP2 interaction. We found that the interaction occurs largely between the PIP2 head group and the cytoplasmic tail domain (CTD) of HA. Mutations of the CTD were made to alter charge (HARE, HARREQ), palmitoylation sites (HAMAY), or a combination thereof (HAREMAY, RREQMAY). MD showed that HARREQ and RREQMAY had the strongest effect on HA-PIP2 interactions through a depletion in the radial distribution function of PIP2 around HA at distances [≤]2.5 nm. FPALM revealed that HA cluster density at the PM was significantly reduced by CTD mutations, with the largest reduction occurring in mutants where the CTD charge and acylation were both altered (HAREMAY). HAREMAY clusters were also found to have larger circularities and perimeters, implying a structural change to the clusters. Mutations in the HA transmembrane domain also caused modest changes to the cluster properties of HA and its co-clustering with PIP2. FPALM showed PIP2 co-clustering with HA was also affected by HA mutations with more free PIP2 localized under HAREMAY clusters. A chemical model of simultaneous HA-PIP2 and PH-PIP2 binding enables interpretation of HA-PIP2 interactions and reveals quantitative differences between PIP2 binding by HA CTD mutants. We conclude that the mechanism of HA-PIP2 interaction consists of at least electrostatic and hydrophobic components. Our insights into the mechanism of HA-PIP2 interaction, and the prevalence of putative PIP2-interacting domains in a number of viral spike proteins suggest it may be fruitful to identify methods of disrupting interactions between phosphoinositides and viral proteins.

In the active loop extrusion model, the cohesin protein complex creates chromatin loops in eukaryotic cells. Extrusion maintains topologically associated domains (TADs), which are contiguous segments of chromatin that preferentially colocalize in space and are typically bounded by CTCF proteins that pause cohesin translocation. Here, we model active loop extrusion with hybrid molecular dynamics - Monte Carlo simulations in entangled flexible linear polymer melts. Intra-chain contact probabilities of polymers with active loop extrusion are enhanced compared to their equilibrium, passive counterparts. Extrusion causes the size of chain segments to be much smaller than in passive melts. While the overlap parameter in passive melts without extrusion monotonically increases with segment length, it is nonmonotonic in active melts and on the order of unity within the parameters of this study. Active loop extrusion suppresses contacts between TADs in favor of intra-TAD contacts. Reduction of overlaps between chain segments dilutes entanglements in active melts. Depending on parameters, active extrusion without TADs may induce more compact conformations than with TADs, due in part to fractal loopy globule-like dynamics. This work suggests that active loop extrusion reduces overlaps between TADs, contributing to effective gene regulation by cis-regulatory elements.

Inefficient endosomal escape is a crucial barrier to intracellular delivery of nucleic acid therapeutics using lipid nanoparticles (LNPs). The use of ionizable lipids (ILs) has significantly improved cargo delivery efficiency, yet the physical basis of their interaction with endosomal membranes and their role in endosomal escape remain unclear. It has been suggested that, as ILs become cationic during endosomal acidification, electrostatic affinity promotes fusion of the LNPs with the endosome. In this paper, we propose an additional mechanism in which ILs are redistributed from LNPs to host membranes, modulating the elastic properties and curvature of the membrane, lowering the energetic threshold for endosome disruption. To test this, we quantified the spontaneous curvature of clinically relevant ILs and ATP-binding lipids and measured the membrane mechanics of giant unilamellar vesicles (GUVs) with an endosome-relevant composition at endosome-relevant pH. Small-angle X-ray scattering (SAXS) measurements reveal that the incorporation of ILs and ATP lipids into endosome-mimetic membranes shifts the spontaneous curvature towards more negative values. Micropipette aspiration experiments indicated a decrease in the apparent area compressibility modulus of membranes doped with ILs and ATP lipids. In addition, membranes showed enhanced fluctuation amplitudes and altered relaxation behavior, consistent with membrane perturbations associated with lipid insertion and pH- or ATP-driven destabilization. Under conditions promoting the partitioning of ILs or ATP-binding lipids, we further observed reduced bending rigidity and increased heterogeneity in membrane tension. Together, these results support a model in which ILs (as well as newly developed ATP-binding lipids) partition into endosomal membranes, softening the membrane and generating local curvature frustration that facilitates endosomal disruption during the natural acidification process. By quantitatively linking lipid composition with changes in membrane elasticity and fluctuation dynamics, this work provides a biophysical framework for understanding how lipid redistribution may contribute to endosomal escape and improve delivery efficiency.

Many membrane proteins contain intrinsically disordered regions (IDRs) that play key biological roles by providing structural plasticity, harboring sites for post-translational modifications, and mediating protein clustering and phase separation. Residue-level molecular models parameterized against experimental data have provided insights into how IDR sequence controls conformational properties and phase behavior in soluble proteins. Here, we extend this modeling framework to membrane-associated IDRs. We adapt a coarse-grained lipid model, iSoLFv2, for four phospholipids and combine it with CALVADOS, a residue-level model for IDRs and multi-domain proteins. Protein-lipid cross-interaction parameters are calibrated to reproduce predicted insertion and orientation of transmembrane proteins with diverse architectures. We validate the resulting model against Wimley-White free energies of transfer of hydrophobic peptides and against an experimentally refined conformational ensemble of a flexible membrane receptor. Finally, we show the applicability of the model to a membrane-associated assembly of signaling proteins. The model provides a computationally efficient framework for studying conformational ensembles and assembly of proteins at bilayer-water interfaces.

Recent experiments show that protein condensates sitting on opposite surfaces of a flat lipid membrane move together and prefer to overlap, even though they cannot touch each other. This points to an indirect, membrane-mediated interaction. Two mechanisms could be responsible: a curvature-induced interaction, which is energetic in origin, and a fluctuation-induced interaction, which is entropic. Here we study both with coarse-grained molecular dynamics simulations, using Cookes implicit-solvent lipid model together with a generic bead-spring polymer model for the condensate. We compute the potential of mean force between two condensates across the membrane. For condensates of the same size, full overlap is unfavorable, and the pair instead settles into a partially overlapping state that bends the membrane into an S-like shape. When the two condensates differ strongly in size, full overlap becomes favorable. We explain this with a simple geometric picture. The condensate wets the membrane as a thin film and imposes curvature only along its rim, while membrane tension flattens the membrane under its interior. The resulting ring of curvature can trap a smaller condensate on the opposite side. We also compare the bending undulations and the effective bending modulus of a bare membrane, a membrane with one condensate, and a membrane with condensates on both sides. A wetting condensate suppresses the undulation modes and stiffens the membrane, but whether this makes overlap entropically favorable remains inconclusive. Our results indicate that the coupling is driven mainly by curvature, and that it depends on the wetting mechanism and on the membrane tension.

A classical study found no excitation transfer when isolated cardiomyocytes were placed side by side, whereas a recent paper reported action potential transfer in carrdiomyocytes placed end to end. We use nanoscale numerical simulations based on the full Poisson-Nernst-Planck equations to investigate whether these apparently opposing observations can be explained by the different geometrical configurations. The computations show that in the end-to-end configuration, ephaptic coupling occurs when the intercellular cleft is sufficiently narrow and a sufficiently large fraction of the sodium channels is localized at the intercalated disc. Coupling is strengthened when the sodium channels are concentrated in fewer clusters and when ionic diffusion within the cleft is reduced. Under these conditions, excitation transfer occurs on a timescale consistent with rapid cell-to-cell activation. Conduction depends biphasically on cleft width and terminates abruptly beyond a critical width. Localization of potassium channels at the intercalated disc has only a moderate effect, whereas gap junctions substantially improve conduction and reduce the relative contribution of ephaptic coupling. In the side-by-side configuration, excitation transfer does not occur under physiological conditions and requires highly flattened cells, minimal separation, and unrealistically strong sodium-channel clustering. The different outcomes of the side-by-side and end-to-end experiments can therefore be explained by the fundamentally different geometrical conditions for ephaptic coupling.

Cellular systems must act robustly to maintain organismal health, including maintaining biophysical properties that allow for appropriate cellular function, and adapting these properties through changes such as those that occur during cell division. However, we still lack tools to measure many of these physical properties with precision in the living cell. For example, the mechanical properties of the nucleoplasm, the fluid-like substance that fills the nucleus, have not been fully characterized. To investigate these properties, we have turned to the fission yeast Schizosaccharomyces pombe (S. pombe), a well-established, genetically tractable model organism that has been used extensively for studying a variety of cell biophysical processes and structures, including the cytoskeleton and cell division. It is an apt system for studying how the nucleus adapts over the course of the cell cycle, since it undergoes closed mitosis, where the nuclear envelope remains intact during cell division. Studying nucleoplasm properties over the course of closed mitosis may help reveal how nuclear volume, shape, surface area expansion, and chromosome segregation are linked and coordinated. To measure nucleoplasm material properties in S. pombe, we have paired Fluorescence Correlation Spectroscopy (FCS) with a photoswitchable fluorophore, enabling fine control over fluorescent intensity inside live cells. We infer material properties from FCS measurements, while the photoswitchable probe enables confocal imaging in conjunction with these measurements, yielding corresponding information about cellular state and dynamics. Interestingly, we find that nucleoplasm material properties do not vary significantly over the cell cycle. Future studies will use this tool to examine how diverse molecular and genetic perturbations alter nucleoplasmic properties, providing insight into how these properties maintain nuclear function and protect genomic integrity over the cell cycle and during development.

Classical neuronal cable theory relies on quasi-static electric field approximations and neglects magnetic induction, Lorentz force coupling, and transient electromagnetic currents, limiting its ability to fully characterize action potential propagation within geometrically branched axons and dendrites. This work develops a coupled Maxwell-electromagnetic cable framework by integrating finite-difference time-domain (FDTD) solutions of Maxwells equations with extended Hodgkin-Huxley and Fitzhugh-Nagumo membrane dynamics, incorporating magnetic gating perturbations, electromagnetic trans-membrane currents IEM, and nanoscale quantum corrections for thin neural segments. Controlled propagation experiments are designed to quantify deviations from standard cable predictions across asymmetric and symmetric axonal bifurcation geometries. Numerical results demonstrate that inductive magnetic effects lower the critical branch radius for junction conduction failure and break symmetric action potential invasion in geometrically identical child branches under external transverse magnetic fields. An electromagnetic corrected geometric ratio GREM is proposed to revise impedance-matching conditions at branch points, accounting for size-dependent axial current imbalance induced by magnetic and displacement currents. Parent axon conduction velocity deviates substantially from the canonical [Formula] scaling law when electromagnetic feedback and quantum charge distributions are included, triggering early signal blockage at large cable diameters. Collectively, this study establishes that quasi-static cable models underestimate electromagnetic corrections to propagation speed, waveform shape, and bifurcation transmission fidelity; the coupled Maxwell-cable framework provides a comprehensive multi-physics tool for modeling electrodynamic signal behavior in complex neuronal architectures.

Maternal thermal strain is associated with adverse pregnancy outcomes, yet fetal temperatures cannot currently be directly measured, limiting quantification of fetal thermal strain. Here, we develop two steady-state models for estimating internal temperatures in a near-term fetus. First, we improve the only previously published human fetal thermoregulation model, deriving a closed-form solution within its simplified uniform-cylinder representation. Second, we introduce a multilayer, anatomically segmented model that resolves tissue-specific temperatures. Both couple the fetal body to central blood pool and amniotic fluid compartments and incorporate a new placenta-umbilical cord heat-exchanger representation. Predictions agree with available intrauterine scalp measurements, with fetal core and head-center temperatures approximately 0.5{degrees}C and 0.8{degrees}C above maternal core, respectively. Physiologically plausible changes in umbilical cord heat-exchanger effectiveness or blood flow increased fetal temperatures by approximately 0.3{degrees}C. These models enable estimation of otherwise inaccessible temperatures, while the multilayer formulation lays a foundation for transient, coupled maternal-fetal thermoregulation modeling.

Reflectin proteins drive dynamic structural coloration in cephalopods by organizing into dense intracellular lamellar structures that dictate local refractive index. While reconstituted reflectins readily undergo liquid-liquid phase separation in vitro, these assemblies frequently undergo dynamic arrest, vitrifying into non-dynamic condensates. Here, we investigate the primary sequence features, post-translational modifications, and heterotypic interactions that regulate the material properties of reflectin condensates within the crowded cellular environment of mammalian HeLa cells. Using confocal microscopy and fluorescence recovery after photobleaching (FRAP), we demonstrate that canonical block copolymeric A-type reflectins readily form dynamically arrested condensates, with the linker blocks primarily responsible for the observed arrest. In contrast, non-canonical B/C reflectin variants exhibit significantly greater fluidity and rapid recovery kinetics. We show that phosphomimetic substitutions progressively fluidize some reflectin condensates. Lastly, we find that heterotypic condensates composed of canonical and non-canonical reflectins in combinations associated with reversible iridescence in squid substantially enhance canonical mobility. Our findings establish a biophysical framework in which phosphorylation and heterotypic mixing cooperatively suppress dynamic arrest, enabling the reversible material transitions required for active cephalopod camouflage and communication.

Cardiolipin (CL) is a phospholipid found in the inner mitochondrial membrane (IMM) where it increases the efficiency of ATP regeneration. We have investigated the hypothesis that this increase may result in part from CL concentrating H+ at the IMM surface through electrostatic interactions as the CL polar head is a dianion at physiological pH. To this aim, we compared the concentrations and movements of H+ at the surface of giant planar phosphatidylcholine (PC) membranes and 20% CL enriched PC membranes by recording their surface pH with the membrane-grafted pH probe fluorescein DHPE. CL enrichment of the membranes increased their surface H+ activity by a ~4 factor. Moreover, we observed non-gaussian spatial H+ concentration profiles with distance from a point H+ source with both PC and CL membranes suggesting that both lipids also induce interactions between probe molecules. A whole bath pH variation revealed that these interactions allow the traveling of reversible acidification fronts with constant speed over the membrane between high and low pH states. A reaction-diffusion model of these observations suggests that membranes support these fronts through a mechanism of autocatalytic (de)protonation of the membrane surface. In mitochondria, these fronts would result in transitions between high and low pH states, the low one having a larger H+ concentration in CL-enriched regions of the IMM. Such an increase at the inner leaflet of the IMM may increase efficiency of the respiratory chain whereas the increase at the outer leaflet may boost the ATP synthase rate.

Cholesterol is a key component of eukaryotic cell membranes, promoting membrane stability and modulating the function of many membrane proteins, including ion channels. In our previous work using purified human voltage-gated proton channel proteins, we showed that cholesterol inhibits the hHv1 channel by altering the conformational dynamics of its S4 segment, the key element that senses membrane voltage to control proton permeation. In the present work, we examined the effects of cholesterol analogs and potential sites in the hHv1 channel mediating cholesterol inhibition using site-directed mutagenesis and docking simulations. Our results showed that desmosterol, the immediate precursor of cholesterol, markedly attenuates cholesterol inhibition. Using single-molecule Fluorescence Resonance Energy Transfer (smFRET), we showed that desmosterol attenuates cholesterol inhibition by promoting the intermediate and open state conformations of the S4 segment. Moreover, we identified multiple residues in the hHv1 channel that are critical for cholesterol inhibition, including Y141A in the S2 segment, which reduces cholesterol inhibition by nearly 3-fold. Our smFRET results showed that the Y141A mutation promotes the intermediate conformation in the S4 segment, which underlies the attenuation of cholesterol inhibition. Consistently, docking simulations also revealed multiple residues spanning the transmembrane domain, rather than clustered within a single localized pocket. Our work identified the key molecular determinant in the hHv1 channel that mediates cholesterol inhibition and also provided a mechanism linking the conversion between demosterol and cholesterol by DHCR24 to pH homeostasis in many cells, such as phagocytes, cardiomyocytes, neurons and microglial cells.

G protein-coupled receptors are major pharmacological targets, yet achieving subtype selectivity remains challenging when closely related receptors share highly conserved orthosteric binding sites. Here, we investigate the molecular determinants governing ligand recognition and unbinding at the adenosine A2A and A2B receptors, two closely related class A GPCRs with markedly different pharmacological profiles. We combine Funnel Metadynamics and adaptive infrequent metadynamics to characterize the thermodynamics and kinetics of three representative ligands: the non-selective antagonist theophylline (TEP), the A2A-selective inverse agonist ZM-241385 (ZMA), and the non-selective full agonist NECA. Across six ligand-receptor complexes, our simulations reproduce experimentally resolved binding modes, predict the unresolved binding poses of TEP and ZMA at A2B, and provide binding free energies consistent with experimental trends. Kinetic simulations further resolve ligand-specific unbinding pathways, metastable intermediates, residence times, rate-determining transitions, and their associated transition-state configurations. Comparison of A2A and A2B reveals how subtle differences within and around their highly conserved orthosteric sites are amplified into distinct thermodynamic and kinetic behaviors. In particular, we identify three major receptor-specific features: differences in hydration and polarity near TM1/TM2/TM7, differences in steric packing and pocket volume at the TM3/TM5/TM6 floor, and a more dynamic network of charged extracellular residues and lipids in A2B that modulates ligand egress. These features rationalize ligand-dependent differences in affinity, residence time, and subtype selectivity, including the preferential stabilization of ZMA-like antagonists at A2A. Overall, our results provide a dynamic atomistic map of the A2A and A2B orthosteric regions and demonstrate how thermodynamic and kinetic information can reveal pharmacologically relevant differences that are not apparent from static structures alone. This framework may support the rational design and repurposing of subtype-selective adenosine receptor ligands.

The mechanical properties of the nucleus are critical for maintaining nuclear integrity and function. We previously showed that chromatin dominates short-extension mechanics whereas lamins provide long-extension strain stiffening. To distinguish the roles of lamin isoforms, micromanipulation nucleus force measurements were performed on isolated nuclei from lamin A/C (Lmna-/-) and lamin B1 (Lmnb1-/-) knockout mouse embryonic fibroblast cells. Lamin A/C knockout does not alter short-extension nuclear stiffness but is essential for strain stiffening at longer extensions. Oppositely, lamin B1 loss reduced short-extension stiffness due to facultative heterochromatin loss while long-extension strain stiffening was slightly increased. Loss of lamin A/C and B1 resulted in similar lamin-chromatin linkers effects as LBR did not change and LAP2{beta} decreased in both. A simulation model of a polymeric lamin shell with stiff lamin A/C and softer lamin B1 subunits can qualitatively recapitulate experimental measurements of lamin knockout cells. Lamin A/C knockout resulted in abnormal nuclear shape but not nuclear blebbing or rupture whereas lamin B1 knockout, similar to other perturbations that cause heterochromatin loss, resulted in increased nuclear blebbing and rupture. This work illuminates the distinct mechanical roles of lamin A/C and B1 in determining nuclear structure and integrity.

The deformability of the red blood cell (RBC) is essential for microcirculatory flow and is profoundly altered in hemoglobinopathies and during infection with Plasmodium falciparum. While many mechanical tests have been developed to probe RBC-mechanics, the dynamics of cell shape recovery following large deformations remains poorly characterized. Here, we integrate microfluidic constriction assays, ultrafast imaging, and computer simulations to quantify time-resolved shape recovery of individual erythrocytes. We show that parasite infection is the primary determinant of RBC viscoelastic behavior. In wild-type (HbAA) erythrocytes, the relaxation time increases progressively from ring to trophozoite to schizont stages, consistent with parasite-induced membrane stiffening and enhanced membrane-cytoskeleton coupling. In contrast, sickle trait (HbAS) erythrocytes exhibit a distinct response: although deformation becomes increasingly irreversible during parasite maturation, the relaxation time after constriction remains largely unchanged. Analysis of a mutant parasite line with enlarged and sparsely distributed knobs revealed a significant increase in relaxation time, demonstrating that knob architecture modulates recovery kinetics. Together, these findings suggest that the coupling between membrane and cytoskeleton, which is strongly changed by the establishment of the knobs during an infection with Plasmodium falciparum, should have a strong detrimental effect on microcirculatory flow, which is however weakened by the sickle cell trait.

The cluster of differentiation 36 (CD36) is a membrane protein with broad physiological roles in health and disease, and its function is regulated in part by phosphorylation. Experimental evidence shows that phosphorylation of Thr92 reduces CD36 affinity for thrombospondin-1 (TSP-1), binding of which initiates antiangiogenic signaling, whereas phosphorylation of Ser237 decreases CD36-mediated fatty acid uptake, with implications for energy metabolism. However, the only available crystal structure of CD36 lacks phosphorylation, and the molecular mechanisms by which phosphorylation regulates CD36 function remain largely unknown. This study provides an atomically detailed computational characterization of CD36 in unphosphorylated and dual phosphorylated states, using molecular dynamics simulations with a total sampling time of 30 microseconds in combination with Markov state models. We present, to our knowledge, the first evidence of a cryptic pocket on CD36 surface that is formed by phosphorylation. This cryptic surface pocket and a loop spanning residues 121-131 form a high affinity binding site for TSP-1 derived ligands, shifting their binding away from the canonical site. We propose that this altered binding provides a molecular basis for the disruption of antiangiogenic signaling upon CD36 phosphorylation. Additionally, our data indicate that, phosphorylation increases helicity and compaction within the helix-loop region spanning residues 296-331, narrowing one of the entrances to the internal cavity and reducing its overall volume. These conformational changes provide a potential mechanistic explanation for the decrease in fatty acid uptake upon CD36 phosphorylation. Our findings provide structural insights that may inform the future design of CD36 modulators and emphasize the importance of targeting phosphorylation induced CD36 conformations in angiogenic and metabolic diseases.

Proteins dynamically switch between a continuum of interconverting conformational states, and understanding these structural dynamics is important for understanding protein function and for developing therapeutics. Molecular dynamics (MD) simulations can provide insight into protein conformational ensembles, but sampling rare conformational states can require substantial computational resources. The recent development of AI-based approaches for generating protein conformational ensembles, such as the Biomolecular Emulator (BioEmu), offers a potential alternative, although it remains unclear whether these approaches can accurately capture the conformational landscapes, especially for special cases such as membrane proteins. Here, we assess the ability of BioEmu to model the conformational dynamics of a model membrane protein, the bacterial rhomboid intramembrane proteases GlpG. We find that BioEmu generates a range of conformations corresponding to both open and closed states of the rhomboid lateral gate, including states associated with different stages of the catalytic cycle. These conformations broadly correspond to states sampled during microsecond-timescale MD simulations, although BioEmu does not reproduce the full conformational landscape observed using MD. BioEmu also samples substantial conformational heterogeneity within the soluble domains of rhomboids, which are highly flexible and poorly represented in experimental structures. Overall, our findings demonstrate that BioEmu can generate plausible conformational ensembles for relatively large, six-and seven-pass membrane proteins, sampling rare states at a fraction of the computational cost of conventional MD simulations. These results suggest that AI-based ensemble generation could provide an accessible approach for exploring membrane protein dynamics and complement conventional molecular modelling approaches.

Focal adhesions maintain force-bearing attachment between cells and the extracellular matrix but can also undergo dynamic remodeling. Their assembly and actomyosin tension are coupled through mechanochemical feedback. The processes underlying this feedback are not instantaneous and therefore involve a time delay. However, how this delayed feedback gives rise to stable adhesion maintenance or dynamic remodeling remains unclear. Here, paired time-lapse measurements of vinculin fluorescence and traction stress revealed distinct local adhesion-force dynamics, including low-fluctuation and recurrent fluctuation patterns. To examine how these patterns could arise, we formulated a minimal mechanochemical model coupling focal adhesion assembly and actomyosin force through delayed reciprocal feedback. The model exhibited stable and oscillatory modes depending on feedback strength, the balance of opposing feedback effects, and the effective feedback delay. Bistability and hysteretic switching also occurred in a subset of parameter space, and the oscillation period followed a power-law relation with the delay. These results suggest that stable adhesion maintenance and dynamic remodeling can emerge from a common mechanochemical feedback architecture.